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ABSTRACT 

We consider trends resulting from two formation mechanisms for short-period super-Earths: planet- 
planet scattering and migration. We model scenarios where these planets originate near the snow 
line in "cold finger" circumstellar disks. Low-mass planet-planet scattering excites planets to low 
periastron orbits only for lower mass stars. With long circularisation times, these planets reside on 
long-period eccentric orbits. Closer formation regions mean planets that reach short-period orbits by 
migration are most common around low-mass stars. Above ~1 Mq, planets massive enough to migrate 
to close-in orbits before the gas disk dissipates are above the critical mass for gas giant formation. 
Thus, there is an upper stellar mass limit for short-period super-Earths that form by migration. If 
disk masses are distributed as a power law, planet frequency increases with metallicity because most 
disks have low masses. For disk masses distributed around a relatively high mass, planet frequency 
decreases with increasing metallicity. As icy planets migrate, they shepherd interior objects toward 
the star, which grow to '^1 M^. In contrast to icy migrators, surviving shepherded planets are rocky. 
Upon reaching short-period orbits, planets are subject to evaporation processes. The closest planets 
may be reduced to rocky or icy cores. Low-mass stars have lower EUV luminosities, so the level of 
evaporation decreases with decreasing stellar mass. 
Subject headings: planetary systems: formation — planetary systems: protoplanetary disks 



\~' ' 1. INTRODUCTION 

> ■ 
\^ ' With nearly 300 known extra-Solar planets, there are 

^\ , now several clear correlations between the properties of 
Cn| ■ the planets and their host stars. The most well known 
CNJ ' trend is the increase in gas giant frequency with host star 
_jj , metallicity (e.g. Fischer & Valenti 2005). Recent radial 
f— ^ ■ velocity surveys suggest that giant planet frequency also 
QQ ■ increases with stellar mass (Johnson et al. 2007). 
f^ , These trends provide tests of planet formation theo- 
?• ■ ries. In the core accretion model for example, gas giant 
. 1^ ' planets form by coagulation of small planetesimals near 
k> , the "snow line" that separates rocky and icy regions in 
; I ' a circumstellar disk. Once icy protoplanets reach a crit- 
d ' ical core mass, they accrete gas rapidly (Pollack et al. 
1996). Cores benefit from extra planet building mate- 
rial provided by enhanced metallicities and an increase 
in disk masses with stellar mass. The model is thus con- 
sistent with current observations (Ida & Lin 2004, 2005; 
Kennedy & Kenyon 2008). 

Gravitational instability (GI) is an alternative forma- 
tion mechanism for gas giant planets, where a relatively 
massive disk cools enough to fragment into Jupiter-mass 
clumps. Although GI operates over a wide range of stel- 
lar masses, there is still debate about predicted trends 
with metallicity (Durisen et al. 2007). Given observa- 
tional biases in the current sample of extra-Solar plan- 
ets, GI cannot be ruled out as a formation mechanism 
(Durisen et al. 2007). 

Core accretion and GI models suggest that short- 
period "hot Jupiters" reside too close to their par- 
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ent stars to have formed in situ. Thus, these planets 
must migrate or scatter from more distant formation 
regions to arrive at their final orbits (Lin et al. 1996; 
Rasio & Ford 1996). A combination of these two mech- 
anisms probably operates to produce the observed dis- 
tribution of extra-Solar giant planets. Scattering can 
reproduce most of the observed eccentricity distribu- 
tion, but has trouble accounting for planets in circu- 
lar orbits at distances too far from their host stars 
for tidal circularisation (Ford & Rasio 2007). Migra- 
tion theories can explain systems with planets in mean- 
motion resonances (Lee & Pealc 2002), but they may 
not reproduce the observed eccentricity distribution (e.g. 
Tremaine & Zakamska 2004). 

With the discovery of the first super-Earths in rela- 
tively short period orbits, migration and scattering re- 
main possible mechanisms for planets to reach these radii 
(Brunini & Cionco 2005; Terquem & Papaloizou 2007; 
Raymond et al. 2007). However, the discovery of low- 
mass planets in systems already harbouring giant planets 
suggests new formation mechanisms (Zhou et al. 2005). 
Because these models require gas giants, they predict 
trends with metallicity and stellar mass for low-mass 
planets similar to those for giant planets. Though some 
low-mass planets may have formed with help from gi- 
ant planets, a flatter metallicity distribution (Udry et al. 
2007) and the absence of giant planets in some low-mass 
planet systems (e.g. G1581 & GJ674, Bonfils et al. 2005, 
2007) indicate other formation mechanisms. 

Here, we consider trends that may arise in forming 
short-period and/or transiting icy/rocky planets in sys- 
tems with no gas giants, over a range of stellar masses. 
The close-in planets that form are therefore the most 
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massive in the planetary system. We first cover some 
background in §2. In §3 we use n-body simulations to 
show that 10 M® planet-planet scattering is unlikely to 
result in transiting planets for all but the lowest mass 
stars. With long circularisation timescales, planets in 
these systems are hard to detect. We consider migra- 
tion scenarios using analytic, semi- analytic and n-body 
models in §4. With migration, short-period low-mass 
planets most likely form around low-mass stars. Above a 
certain stellar mass, it is hard to form any short-period 
planets without giant atmospheres. Trends with nietal- 
licity depend on the disk mass distribution. Migration 
to short-period orbits results in significant amounts of 
material being shepherded inward, which affects the fi- 
nal structure of these systems. We discuss our results, 
subsequent planetary evolution, and conclude in §5. 

2. BACKGROUND 

2.1. General Picture 

Planets form in circumstellar disks. Therefore disk 
structure plays a key role in setting the final configu- 
ration of planetary systems. In most planet formation 
models, disk structure is characterised by an outwardly 
decreasing radial surface density profile. This profile usu- 
ally includes an increase in surface density at the "snow 
line," where the temperature becomes low enough for 
water to freeze. 

Planets form by accumulating solids in the disk. 
Therefore the expected increase in surface density at the 
snow line is often associated with the formation of gas 
giants like Jupiter. Forming Jupiter requires the rela- 
tively rapid growth of a ^^5-10 M,^ icy core, followed by 
a period of gas accretion (Pollack et al. 1996). Gas ac- 
cretion must be complete before the gas disk disperses 
in '--'3Myr (e.g. Haisch et al. 2001). In the minimum 
mass Solar nebula model (MMSN, Weidenschilling 1977; 
Hayashi 1981), forming the icy core rapidly requires 
factor of 5-10 surface density enhancements relative to 
the terrestrial region (Lissauer 1987; Pollack et al. 1996; 
Thommes et al. 2003) . This factor is larger than the fac- 
tor of 2-3 enhancements expected from Solar abundances 
(Asphmd et al. 2005), or suggested by comet composi- 
tion (Kiippers et al. 2005) , and the factor of ^4 derived 
in the original MMSN model (Hayashi 1981). 

The need for larger surface density enhancements 
inspired "cold finger" disk models, which produce 
much larger snow line enhancements in a relatively 
narrow (^AU) radial region near the snow line 
(Stevenson & Lunine 1988; Cuzzi & Zahnle 2004). In 
this picture, a circumstellar disk has an initial equilib- 
rium state with the water vapour (ice) concentration 
decreasing (increasing) beyond the snow line. As the 
disk diffuses and advects, water continually condenses 
from gas passing beyond the snow line, thus enhancing 
the local surface density of solids, and removing vapour 
phase water from the inner disk. Sublimation of plan- 
etesimals that drift inside the condensation radius by gas 
drag enhances this effect: the surface density beyond the 
snow line increases when water vapour from the subli- 
mated planetesimals diffuses back outside the snow line 
(Cuzzi & Zahnle 2004). 

The first cold finger models predict a factor of 
~10-100 increase in the surface density of icy mate- 
rial in a relatively narrow region near the snow line 



(Stevenson & Lunine 1988; Cuzzi & Zahnle 2004). Us- 
ing a more complex global disk model, Ciesla & Cuzzi 
(2006) suggest surface density enhancements closer to 10 
than 100. In their simulations, the enhancement regions 
are several AU wide at half the maximum planetesimal 
surface density. 

The main differences expected for planet formation 
models using cold finger instead of MMSN disks are 
threefold. Due to the nature of the surface density en- 
hancement: (i) fewer large planets form, (ii) large planets 
form in relatively low-mass disks, and (iii) planets form 
from material with much higher ice/rock ratios. In ad- 
dition, material lost to inward planetesimal drift by gas 
drag (Thommes et al. 2003) may be returned to the cold 
finger region, allowing continued growth. Reducing the 
removal of drifting planetesimals enhances growth rates, 
and allows formation of more massive icy planets. 

2.2. Mathematical Formalism 

In the standard coagulation model, planets grow in a 
circumstellar disk through repeated collisions and merg- 
ers of smaller objects (Safronov 1969). First, roughly 
km size planetesimals form rapidly, whether by coagu- 
lation (e.g. Weidenschilhng 2000) or direct collapse (e.g. 
Goldreich & Ward 1973). Little knowledge of which pro- 
cess dominates means the size distribution of the first 
planetesimals is poorly constrained. Planetesimals ini- 
tially grow through a rapid phase of "runaway" growth 
(Kokubo & Ida 1996). During the period of "oligarchic" 
growth that follows (Kokubo & Ida 1998), protoplanc- 
tary growth rates depend on the surface density of plan- 
etesimals (7s, the local orbital frequency £7, the gravita- 
tional reach of the growing protoplanet, and the random 
velocities of the smaller planetesimals (Inaba et al. 2001) 

A/pi oc o-g Th r2 Pcoi(e, i) ■ (1) 

Here rn = a (A'/pi/3Af*) ' is the Hill radius, and a is 
semi-major axis. The eccentricity e and inclination i 
are in units of the growing protoplanets Hill radius (i.e. 
e = e/y'n). The collision probability Pcoi largely de- 
termines how growth proceeds: growth is fastest when 
planetesimals are small enough (<lkm) to be damped 
by gas drag (e.g. Rafikov 2004). In this "shear domi- 
nated" regime when e and i are <1, growth depends on 
Keplerian shear in the disk, rather than objects random 
velocities. Growth slows strongly with increasing radial 
distance, because J7 ex a^'^/^ and CTs oc a^^ , where 5 ~ 1- 
1.5. 

Eventually, protoplanets accrete most of the nearby 
material and reach the "isolation" mass (Lissauer 1987) 

Mi.„ = ^'"" '° .i . (2) 
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Numerical simulations indicate that isolated oligarchs are 
spaced at 2i?i?H ^ 8i?H intervals (e.g. Kokubo & Ida 
1998). In the terrestrial region around the Sun, the iso- 
lation mass is ~0.1 M®, and the timescale for Earth for- 
mation by the chaotic growth that follows is ^10-100 Myr 
(e.g. Kenyon & Bromley 2006). 

Further out in the disk, larger isolation masses allow 
formation of gas giant planets. The critical core mass for 
gas accretion depends on opacity and planetesimal accre- 
tion rates, but is >10 Af^ (e.g. Ikoma et al. 2000; Rafikov 
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2006). This mass is reached more easily further out in 
the disk because Miso increases with a. However, growth 
slows rapidly with increasing radial distance; thus, there 
is an optimum region where cores are massive enough 
to accrete gas and to form giant planets before the gas 
disk is dissipated (Kennedy & Kenyon 2008). This re- 
gion is sufficiently far from the star that in situ formation 
of "hot-Jupiters" is unlikely, thus motivating theories of 
migration and scattering. 

2.3. Migration 

Type I migration is a potential barrier to the 
formation of both terrestrial and giant planets 
(Goldrcich & Tremaine 1980; Ward 1997; Tanaka et al. 
2002; Papaloizou et al. 2007). When protoplanets reach 
near an Earth mass, the excitation of spiral density waves 
in the gaseous disk causes planets to experience a torque, 
and migrate inward. The timescale for a planet to spiral 
into the central star is (Tanaka et al. 2002) 



(2.7+1.16) 



_i {M^MqY K" 



(3) 



where h w 0.05 is the disk aspect ratio, and the stellar 
mass Af* is in units of Solar masses. For a planet of 
mass Mpi = 1 Mq in a disk with Cgas — 1700 g cm~^ at 
1 AU around a Solar-mass star, Tmig = 1.6 x 10^ yr. Be- 
cause this timescale is shorter than the ^^3 Myr disk life- 
time (Haisch et al. 2001), and comparable with growth 
timescales, type I migration theory conflicts with terres- 
trial and giant planet formation in the Solar System (but 
see Chambers 2006b). 

Recent work suggests a reduced migration efficiency 
can resolve this problem (Ida & Lin 2008). This "off- 
set" applies to planets <15M© (D'Angelo et al. 2002, 
2003) and arises from corotation torques by coorbital 
material (Masset et al. 2006a). Other ways of reduc- 
ing (and even reversing) type I migration rates include 
turbulence arising from the magnetorotational instability 
(e.g. Nelson k. Papaloizou 2004), and eccentricity driven 
by planet-planet interactions (Papaloizou & Larwood 
2000). 

If planets do not fall onto the central star, migration 
is a possible mechanism for producing planets on short- 
period orbits (Lin et al. 1996; Brunini & Cionco 2005; 
Terquem & Papaloizou 2007). 

2.4. Scattering 

Planet-planet scattering can also produce planets with 
short-period, or low periastron [q) orbits. Originally pro- 
posed to explain hot-Jupiters (Rasio & Ford 1996), this 
scenario has not been applied to low-mass planets. 

Scattering favours giant planets on short-period orbits. 
When a gas giant scatters into a low periastron orbit, 
tidal interaction with the star can circularise the orbit 
on reasonable timescales, with a ~ 2q (Rasio & Ford 
1996). For lower mass planets, long circularisation 
timescales make circular orbits unlikely (Raymond et al. 
2007). However, if the initial scattering region is suf- 
ficiently close, as for low-mass stars, detection of low- 
periastron eccentric planets is possible. 

We now consider two different scenarios that form 
short-period and/or transiting low-mass planets that be- 
gin growth near the snow line, across a range of stel- 
lar masses. When the snow line enhancement is small. 



many planets migrate toward close orbits. This sce- 
nario has already been studied for Solar-mass stars by 
Terquem & Papaloizou (2007). Here, we instead con- 
sider cold finger type disks, where a few planets forming 
near the snow line dominate others forming elsewhere in 
the disk. We first consider a scattering scenario result- 
ing from in situ growth, and then a migration scenario. 
We defer discussion of subsequent planetary evolution in 
final orbits to §5. 

3. SCATTERING 

Planet-planet scattering is a likely outcome of oli- 
garchic growth. In migration scenarios, protoplanets in- 
teract strongly with the gas disk, and they migrate to 
close-in orbits. However, if the gas disk disperses be- 
fore planets have time to migrate, or if migration re- 
sults in no net inward movement, planets form in situ. 
During oligarchic growth, protoplanets grow on orbits 
near the limits of dynamical stability, with damping pro- 
vided by small bodies (e.g. Stewart & Wctherill 1988; 
Kokubo & Ida 1998). At later stages near isolation, their 
orbits can become unstable as remaining small bodies 
are accreted (Goldrcich et al. 2004; Kenyon & Bromley 
2006). 

When planets start interacting dynamically, the 
boundary in semi-major axis between stable and unsta- 
ble configurations is very sharp. Thus, two planets with 
orbits that become too close experience the sudden on- 
set of a dynamical instability caused by close encounters 
(Gladman 1993). 

In previous studies of giant planet scattering, plan- 
ets begin at ~AU distances from the central star, with 
spacings just inside the stability limit. After many inter- 
actions, one planet sometimes attains a highly eccentric 
orbit with a small periastron distance (e.g. Rasio & Ford 
1996; Ford & Rasio 2007). Tidal interaction with the 
central star then circularises the orbit with a ~ 2q. 

While tidal forces can circularise gas giant orbits, the 
timescales for 1-10 Mq planets on highly eccentric orbits 
are long (>Gyr, Raymond et al. 2007). Although these 
planets maintain eccentric long-period orbits, transits are 
possible in favourable circumstances. Because planets 
form at shorter orbital periods around low-mass stars, 
these provide the best opportunity for transit observa- 
tions. 

Cold finger disks provide an ideal environment for oli- 
garchic growth followed by planet-planet scattering. The 
width of the cold finger region allows several protoplanets 
to form (Ciesla & Cuzzi 2006). Once protoplanets reach 
isolation, further chaotic growth may occur if their es- 
cape velocity Vosc is less than the local Keplcrian velocity 
wk ("^ = t'esc/wK, Goldrcich et al. 2004). In the terres- 
trial region of Solar-type stars, TZ ~ 1/4. For gas giants, 
7^ > 1. For A/pi = 10 Me with density p = 4.5 g cm"^, 
TZ K, 1.3 outside the snow line. Thus, '^10 Mq protoplan- 
ets present an approximate division between coalescence 
and scattering/ejection, and an order of magnitude esti- 
mate of the maximum planet mass. This mass is similar 
to the minimum needed for gas accretion, so scattering 
of super-Earths to close-in orbits appears difficult. 

For less massive stars, scattering to low periastron or- 
bits is easier. At fixed a, smaller wk leads to larger TZ 
and a greater chance of scattering. However, the snow 
line also moves inward as stellar mass decreases (csnow oc 



KENNEDY & KENYON 



TABLE 1 

Scattering simulation outcomes, and fraction with 
low periastra. 



M^Mq) 


collisions 


ejections 


survival 


g<0.1AU 


0.25 


95% 


2.5% 


2.5% 


4% 


0.5 


95% 


0% 


5% 


0% 


1 


94% 


0% 


6% 


0% 


2 


84% 


0% 


16% 


0% 



M^i-2, e.g. Ida & Lin 2008; Kennedy & Kenyon 2008), 
so scattering remains difficult. For Osnow oc Af^,, 
Wif(asnow) is constant for different stellar masses. How- 
ever, for a fixed time period, a greater number of con- 
junctions for low-mass stars allows dynamical evolution 
to greater eccentricities. 

3.1. Scattering Simulations 

To measure the likelihood of planet-planet scattering, 
we performed simulations over a range of stellar masses 
with the MERCURY integrator (Chambers 1999). We 
initialised integrations with two 10 M® planets spaced 
near the Hill stability criterion to ensure close encoun- 
ters (Gladman 1993). This planet mass is an approxi- 
mate maximum mass before cores accrete gas to become 
gas giants, and thus offers the best chance for scattering 
over coalescence. To represent a linearly stellar mass de- 
pendent snow line, the inner planet was placed at ain = 
3 Mj, AU. The outer planet begins at a random a in the 

range 0.9-1 ain (1 -I- Adit), where Adit = 3 (A/pi/A/^) 
(Gladman 1993; Ford & Rasio 2007). Both planets be- 
gin in circular orbits with random inclinations less than 
3°; the remaining orbital elements are chosen randomly. 
Simulations were run with a 5 day timestep for 1 Gyr 
around stars of 0.25, 0.5, 1, and 2Mq, or halted earlier 
in the case of collisions (we assume perfect mergers) or 
ejections. A total of 520 simulations were run, 130 for 
each stellar mass. 

3.2. Scattering Results 

The simulations result in three different outcomes: col- 
lisions, ejections, or survival of both planets for 1 Gyr. 
No planets achieved periastra low enough to fall onto 
the central star. Most (> 85%) simulations resulted in 
collisions (Table 1). Some systems survived for the full 
simulation. The only ejections were for 0.25 Mq. 

With so few systems remaining after 1 Gyr, we use the 
smallest periastron distance reached in each simulation 
to characterise the success of planet-planet scattering, 
shown in Figure 1. As expected, the closer snow line 
distance for the 0.25 M© allows smaller periastra after 
scattering. 

For simulations of 0.25 M© stars, 5/130 (4%) plan- 
ets reach periastra less than 0.1 AU. A shorter orbital 
period allows many more conjunctions. Thus, systems 
evolve further than for more massive stars. For the three 
ejections, the lowest periastra were reached just before 
a series of close encounters, which resulted in the ejec- 
tion. For the three surviving systems, the lowest pe- 
riastra were reached near the end of the integrations. 
These orbits have eccentricities saO.5, and semi-major 
axes «0.5AU, corresponding to an orbital period of 
around 260 days. Circularisation times for these planets 




Fig. 1. — Smallest periastra of scattering simulations for all non- 
stable orbits. Black filled bars: simulations in which both planets 
survived until 1 Gyr. Grey bars: ejections (orbits going beyond 
1000 AU). Unfilled bars: collisions/mergers. The dashed line shows 
the starting a (and q) of the inner planets (circular) orbit. Planets 
are scattered to lower periastra for low-mass stars because of more 
conjunctions, and starting closer to the central star. 



are ~10Gyr (Goldrcich & Sotcr 1966; Raymond et al. 
2007). 

Therefore, in the case of 10 M^ planet-planet scatter- 
ing, only the lowest mass stars have planets with perias- 
tra close enough for transiting orbits. However, long cir- 
cularisation timescales mean these planets will likely re- 
main on highly eccentric orbits, with periods long enough 
to make radial velocity and transit detections difficult. 

4. MIGRATION 

We now turn to a migration scenario. In this picture, 
the largest protoplanets form near the snow line by oli- 
garchic growth. Once they reach masses of ^I-IOM^, 
these icy protoplanets migrate towards the central star. 
In a cold finger disk, migrating icy objects dominate 
smaller interior rocky protoplanets. Long chaotic growth 
timescales mean that as the icy object migrates through 
the terrestrial region, interior rocky objects are not ac- 
creted. They are instead scattered outward or shep- 
herded inward. Shepherding — where interior objects are 
captured into mean-motion resonances — results in rocky 
protoplanets being pushed inward ahead of the migrat- 
ing icy planet. These smaller objects merge to form large 
rocky planets, which are eventually accreted by the larger 
icy migrator or survive on an interior orbit. We assume 
that all objects halt their migration when they reach the 
inner edge of the gas disk, at ~10 stellar radii. 

Our goal is to calculate the growth and migration of 
individual protoplanets in this scenario. If several pro- 
toplanets migrate, whether they do so as a set in res- 
onant orbits (Terquem & Papaloizou 2007), or succes- 
sively (e.g. Daisaka ct al. 2006; Ida & Lin 2008), the re- 
sulting trends are similar in our picture. The main dif- 
ference between outcomes is the number of icy planets 
on short-period orbits. The trends our models predict 
depend largely on disk and planet properties, not multi- 
plicity. 

When an icy protoplanet migrates, it only interacts 
dynamically with interior objects. Because collision 
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cross sections are essentially geometric, the timescale 
for growth is much longer than the migration timescale. 
For example, the migration timescale for an Earth-mass 
planet at 1 AU (several 10^ yr, or several 10^ yr if mi- 
gration is less efficient) is much smaller than the chaotic 
growth timescale (rchaotic ~ pRpi/cTs^ ~ fO^yr, where 
Rpi is the planet radius, Goldreich et al. 2004). Thus, the 
migrating protoplanet does not accrete terrestrial proto- 
plancts; outward scattering or inward shepherding arc 
the most likely outcomes. 

The evolution of interior protoplanets depends on their 
random velocities. Chaotically growing objects with high 
eccentricities are scattered outward by the migrating pro- 
toplanet. These may interact with another migrating 
protoplanet or resume chaotic growth. If interior ob- 
jects have finished chaotic growth and are damped by 
the gas disk onto more circular orbits, shepherding by 
capture onto resonant orbits is possible. Shepherded ob- 
jects merge and form rocky planets as their orbits are 
pushed together by the migrating icy protoplanet. Shep- 
herding by giant planets undergoing type H migration 
has been proposed as a way to form super Earth-mass 
planets (Zhou et al. 2005). However, studies have yet to 
consider shepherding by super-Earths undergoing type I 
migration. 

While some planets are stranded at intermediate radii 
as the gas disk dissipates, most planets that begin to mi- 
grate reach the inner disk edge, and might fall onto the 
star. Because the torque on the migrating planet changes 
when the disk gas surface density profile varies rapidly, as 
happens at the inner disk edge, this fate may be avoided 
(Tanaka et al. 2002). Here, corotation torques affect mi- 
gration, and allow for planets to cease migration before 
reaching the stellar surface (Masset et al. 2006b). In our 
migration simulations, we therefore assume migration 
stops inside the inner disk edge (Terqucm & Papaloizou 
2007). 

In the rest of this section, we consider three models 
that explore different aspects of the migration scenario, 
and observable trends that probe stellar and disk prop- 
erties. We consider the simplest scenario — when growth 
is so fast that planets reach isolation before migration 
begins — with an analytic model in §4.1. As the plan- 
etesimal size increases, growth slows; the timescale be- 
comes comparable to that for migration. The assumption 
made in the analytic model no longer applies, and we 
use a semi-analytic model to study concurrent growth 
and migration in §4.2. Finally, we use n-body simula- 
tions in §4.3 to show the shepherding effects migrating 
super-Earths have on terrestrial material. 

4.1. An Analytic Approach 

If we assume that protoplanets reach isolation before 
migration starts, then we can create a simple analyti- 
cal model for our migration scenario. At isolation, pro- 
toplanets have a known migration timescale, which is 
shorter than the disk lifetime if they are to reach the 
central star. To remain in the super-Earth mass regime, 
the mass of a protoplanet is smaller than the critical 
core mass for gas accretion. Because the isolation mass 
changes with surface density — and thus with disk mass — 
only a certain range of disk masses satisfy these con- 
ditions for fixed stellar mass. To consider a range of 
different stars, we also consider how the snow line — 



where these migrating planets form — changes with stellar 
mass. The range of disk masses that satisfy the condi- 
tions changes with stellar mass, resulting in potentially 
observable trends that test migration models. 

To begin, we adopt a relation for the surface density of 
solid material in the disk. In the standard MMSN model. 



O'Q /ico a, 



-s 

AU ' 



(4) 



where ctq = 8g cni^^, S = 1-1.5, and cau is a, in units 
of AU. The factor /ico ^ 2-3 is the enhancement from 
ice condensation beyond the snow line. This disk has a 
mass --O.OIMq. 

To generalise this relation, wc add terms to account 
for differences in disk mass and metallicity around 
stars with a range of masses. Disks around young 
stars have a large dispersion in mass (Natta et al. 2000; 
Andrews & WiUiams 2005, 2007). Setting the disk mass 
Mdisk oc rjM^ allows us to treat the observed trends with 
stellar mass — Mdisk oc M-f , with (3 « 1 — and a range (77) 
of disk masses at fixed stellar mass. Adopting a factor 
^A oc lO^'^/^l for the metallicity of the stars and the disk 
yields 

as ^ ao V. f ice MM^al^. (5) 

For simplicity, we combine /ice and M. into a single factor 
^ = ficc-M, which quantifies the enhancement of solid 
material relative to gas where these planets form. For a 
cold finger disk, we use /ice = 10. Thus, for typical ranges 
in A^ (^^1/3-3) and /ice (2-10), the plausible range of A 
is 0.6-30. We concentrate on higher A, because these are 
cold finger disks. 

For the surface density of the gas disk, we set dg = 
lOOiJs/A. Thus, the gas mass depends on 77 and S, and is 
independent of metallicity and the enhancement in ices 
at the snow line. We adopt S ~ 3/2. 

How the snow line varies with stellar mass is uncer- 
tain. The existence of gas giant planets suggests that 
the stages of planet formation up to isolation occur while 
the gas disk is still present. During these stages the snow 
line distance is set by viscous accretion of the gas disk. 
If the accretion rate onto the star is M oc M^~^, then 

flsnow oc Mi,' ~ ' (Kennedy & Kenyon 2008). Later, 
when the star has reached the main-sequence and the gas 
disk has been dissipated, the main-sequence luminosity 
is more important and flgnow oc M^ (Ida & Lin 2005). 
Because we model oligarchic growth, and M oc M^ 
(MuzcroUe et al. 2005), we adopt the snow line distance 
Osnow = 2.7Mi,A\J. Variation of the snow line with 
time and stellar mass is a key component of planet for- 
mation models that consider a range of spectral types 
(Kennedy & Kenyon 2008). 

Substituting our adopted surface density into the iso- 
lation mass yields 



Mis 



CTsa" 



\3/2 



(77/iccA^M.ai/2)3/2 



(M.) 



1/2 



m: 



1/2 



(6) 



The isolation mass increases with any parameter that 
increases the surface density. The increasing disk mass 
with stellar mass (Mi, in numerator) is stronger than the 
decreasing Hill radius (M^, in denominator). Thus, at 
fixed a the isolation mass increases with stellar mass. 
For our scenario, we arc interested in planets that form 
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at the snow line, so the changing snow hne distance (a = 
Osnow oc M*) makes the stellar mass dependence stronger. 
Substituting a = 2.1 M^, AU, equation (6) yields the iso- 
lation mass at the snow line for a range of stellar and disk 
masses, and metallicities and snow line enhancements 



M,, 



0.12 {/^rjf'^Ml'^M^ 



(7) 



Applying the same approach to type I migration yields 



1? Mpi n a. 



Ml'' h^ /, 



mig ' 



gas 



MpiTy 



(8) 



where the offset fmig, allows us to consider reduced migra- 
tion rates. At fixed a, migration takes longer as stellar 
mass increases, and speeds up as planet mass increases. 
If planet masses vary less strongly with radial distance 
than Mpi oc a, then the migration timescale increases 
outward, and planets cannot catch up to interior ones. 
Even with isolated objects (Miso oc a'^/^), planets may 
not catch up to interior ones due to the strong slowing of 
growth with semi-major axis. At the snow line distance, 
migration slows even more strongly with increasing stel- 
lar mass due to lower gas density and slower orbital pe- 
riods at larger radii. 

Again substituting a = 2.7Afj,AU, the timescale to 
migrate from the snow line to the star is 



9.1x10^-'''°'^^^* 



3/2 



A/pi V 



■yr, 



(9) 



where we have set h = 0.05 (e.g. Papaloizou et al. 2007). 
At fixed Mi,, massive planets in massive disks migrate 
to the inner disk edge fastest. The migration timescale 
increases with Mj, because the snow line is further away. 
If the migration time is shorter than the disk lifetime 
(i.e. Tniig < Tdisk ^ IMyr), then protoplanets reach 
short-period orbits. This inequality leads to 



Mpi> 



0.91 /.nig M: 



3/2 



Mff 



(10) 



This result yields the minimum mass for a planet to mi- 
grate to a close orbit. Substituting the isolation mass 
(eq. 7) for Mpi and solving for -q gives a lower relative 
disk mass limit of 



ri > mow 



2.2/, 



2/5 



A/y'°A3/5 



(11) 



Disks more massive than this jy form protoplanets mas- 
sive enough to migrate to short-period orbits before the 
gas disk dissipates. Planets in slightly less massive disks 
still migrate, but are stranded at intermediate radii as 
the disk disperses. 

The critical '^lO Af^ core mass for gas accretion pro- 
vides an upper limit for the protoplanet mass. Solving 
A^iso < 10 Mq for r] yields 



V < rjhi = 



30 



18.6 



fs.e 



<30 
>30 



(12) 



where the additional constraint of a reasonable disk mass 
sets Tj < 'iQ (A/disk ^ 0.3 M^,) as an upper limit (e.g. 
Ida & Lin 2005). Because we assume growth is fast, 
planctesimal accretion drops significantly at later stages. 




Fig. 2. — The range of planet masses that reach close-in orbits 
as a function of stellar mass. The thicker grey line line shows re- 
sults from the analytic model using equations (7), (11), and (12). 
Other lines show upper and lower mass limits for a range of plan- 
ctesimal sizes from the semi-analytic model (see §4.2) for A = 10: 
r = 10m (plain), 100m (dashed), 1km (dot-dashed), and 10km 
(dotted). The range of planet masses reaching short-period orbits 
decreases with increasing stellar mass because the snow line dis- 
tance is greater. The inset panel (same axes) shows how different 
snow line relations affect the model (using agnow = 2.7M^AU). 
Lines are for a = 1/2 (dotted), 1 (plain), and 2 (dashed) A more 
strongly varying snow line distance (a = 2) yields smaller Miao 
(due to smaller iJu) S't niuch closer snow line distances as stellar 
mass decreases. 



The core mass for gas accretion is then somewhat smaller 
(Ikoma et al. 2000; Rafikov 2006). 

The two limits on disk mass yield a simple relation 
between the stellar mass, migration offset, and enhance- 
ment factor. Equating ?7iow and r^hi. 



-^^Tt.max — 



7.3 



(/migA)3/« 



(13) 



in units of Solar masses.^ This equation has a sim- 
ple physical interpretation. For massive stars (Af* > 
A/jr.max), the only protoplanets massive enough to mi- 
grate to the central star before the gas disk disperses 
are above the critical core mass for gas accretion. These 
planets therefore become gas giants, rather than forming 
hot super-Earths. For lower stellar masses, the closer 
snow line distance allows planets smaller than the criti- 
cal core mass to migrate to the host star. Thus, A^^^^max 
is the maximum stellar mass for hot super-Earths pro- 
duced by type I migration. 

Making an estimate of Afj^^max requires an assumed 
/mig and A. For Solar metallicity A^ = 1, and a cold 
finger enhancement /ice = 10-20, A = 10-20. For a mi- 
gration offset /niig = 10, Afj^^niax ~ 1 Af . Transit and ra- 
dial velocity surveys routinely probe these stellar masses. 
Independent of the disk mass distribution, this result is 
therefore a simple testable prediction of hot super-Earth 
formation by type I migration. 

Figure 2 shows the range of planet masses that reach 
short-period orbits for a range of stellar masses. For 
the analytic model (thick grey lines) the upper limit is 
constant at 10 Af^. The lower limit decreases as stel- 
lar mass and snow line distance decrease. The expected 

^ Equating (11) and (12) has two solutions for Mi, because of 
the upper limit of 30. The other solution is at A/* far too small to 
be interesting. 
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Fig. 3. — Probability distributions for power law (left), and Gaussian with /J. = 1 (right) disk mass distributions of <1OM0 planets as 
a function of stellar mass. For /mig = 10, lines are for A = 0.5, 1, 2, 5, 10, 20, and 50 from right to left where curves cross the x-axis. 
Thus M*_max decreases with increasing A, and is independent of the disk mass distribution). Each plot is arbitrarily normalised to 1 for 
the most likely A at M* = 0.1 Mq. 



range of planet masses decreases with increasing stellar 
mass, while the average mass increases to 10 M®, where 
the lines meet at M^^^max = 1-3 M©. 

In addition to this maximum stellar mass, we can de- 
rive the probability of forming hot super-Earths around 
stars with M^, < il/^.i„ax- This estimate requires an 
adopted distribution of rj (i.e. disk masses). If rela- 
tive disk masses (Mdisk/M^,) are distributed as a power 
law with index ~— 1.75 (Andrews & Wilhams 2005), the 
(relative) probability of forming a close-in planet as a 
function of stellar mass for a given A is 

Pp(A/„A)(x / 77-1-75^77. (14) 



Via 



Alternatively, 
some "typical'' 

Pg(Af.,A) 



disk masses may be distributed around 
relative disk mass (e.g. Ida & Lin 2005) 

exp I - ' •"^' , • ' dri (15) 



I7hi 



Via 



2s2 



where we choose the standard deviation s = 1. This 
distribution is plausible because opacities may under- 
estimate disk masses by as much as an order of mag- 
nitude, due to mass locked up in boulder size objects 
(Andrews & Williams 2007) . Therefore mm observations 
see disks not only with a range of masses, but in a range 
of evolutionary states. Unlike the case for giant plan- 
ets, there is no observational anchor point, so we present 
these results as relative probabilities. 

For a range of A, the left panel of Figure 3 shows 
the probability distribution for the power law disk mass 
distribution with /mig = 10. Results are similar for 
/mig = 1, with the main difference that M^^^max is higher 
(eq. 13). Higher A are most relevant here because 
low values describe MMSN disks, which result in many 
similar-mass migrating planets originating from a wide 
range of radii. The point where lines break and decrease 
toward lower stellar masses is caused by the maximum 
disk mass condition 77 < 30. In these cases the maximum 
short-period planet mass is not set by gas accretion, and 
is <1OM0. This limit applies when A < 5 for the lower of 
the stellar masses we consider, so does not apply to cold 
finger disks with /ice ^10 unless they have metallicity 
M < 0.5. 



At the lowest stellar masses, there is a clear increase 
in planet frequency with A. With a power-law distribu- 
tion of disk masses, the most common disks are the least 
massive; these require large A to allow them to form 
planets massive enough to migrate (and satisfy condi- 
tion 11). Near M^^^niax, there is an optimum A, which is 
a balance between the likelihood of different disk masses 
and the A needed to form close-in planets from those 
disks. At M*,,nax, the only planet that reaches a short- 
period orbit has A/pi = 10 Mq. Therefore the range of 
short-period planet masses decreases up to Af^, „iax- The 
average planet mass increases with stellar mass. 

The right panel of Figure 3 shows the probability dis- 
tribution for the Gaussian distribution with /i = 1. The 
most common disk mass is thus '^0.1 Ah. As A increases, 
the probability of forming a short-period planet decreases 
once the disk mass distribution is not truncated by the 
condition 77 < 30. In contrast to the power law distri- 
bution, the low-mass disks requiring large A are uncom- 
mon. Thus, as A increases, isolation masses are pushed 
over the gas accretion mass, and the likelihood of form- 
ing close-in <lOAf0 planets decreases. While the curves 
are different from the left panel, the point M^^max is the 
same for a given A. With ^ = (i.e. distributed about 
A^disk = 0.01 Mi,) the probability distribution is qualita- 
tively similar to the power law disk distribution. 

In summary, the simple analytical model yields 
testable predictions for an ensemble of super-Earths that 
migrate into short-period orbits from the snow line. For 
reasonable input parameters, we predict a maximum stel- 
lar mass '--^l Mq for stars with close-in super-Earths. 
If circumstellar disks tend to have similar snow line 
enhancements, this maximum mass decreases with the 
metallicity of the host star. For a range of stellar masses, 
the frequency of hot super-Earths depends on the ini- 
tial distribution of disk masses. For a power-law (Gaus- 
sian) distribution of disk masses, the model predicts more 
(fewer) hot super-Earths around more metal-rich stars. 

To give the these trends some context, the first tran- 
siting low-mass planet orbits a star with sub-Solar mass 
and metallicity (GJ436b, Gillon et al. 2007). The current 
sample of low minimum-mass planets also indicates a 
flatter metallicity distribution than exists for giant extra- 
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Solar planets (Udry ct al. 2007). While both disk mass 
distributions suggest that low stellar mass host is likely, 
the power law distribution argues against a low metallic- 
ity host. The Gaussian disk mass distribution, centered 
on a relatively high disk mass is consistent with an in- 
creasing giant planet frequency with metallicity, and a 
flatter or decreasing frequency for lower mass planets. 

Disks with rj > 77hi form gas giants. Their relative 
probabilities can thus be calculated by integrating equa- 
tions (14) and (15) from ?7hi to 30. However, because 
?7iow only weakly depends on M^,, giant planet frequency 
is roughly some constant minus the hot super-Earth fre- 
quency (i.e. generally increases with M*). This trend is 
essentially the result arrived at by previous theoretical 
studies (e.g. Ida & Lin 2005; Kennedy & Kenyon 2008), 
and is at least qualitatively consistent with the observed 
trend (Johnson et al. 2007). 

In constructing the above model we simplified some pa- 
rameters, and assumed values for others. We now briefly 
consider model sensitivity to these, and whether obser- 
vations may constrain them. The most uncertain sim- 
plification is how the snow line distance varies with stel- 
lar mass. Within our framework, relaxing the distance 
to flsnow ~ 2.7 M" AU results in changes to Equations 
(12), (11), and (13) for a = 1/2-2 (Fig. 2 inset). A 
more strongly varying snow line distance (a = 2) yields 
much closer Osnow and smaller Mjso (due to smaller Rh) 
for low mass stars. A more complex snow line model 
could include how Csnow varies with Mdisk at fixed stel- 
lar mass, or some time dependence (e.g. Ciesla & Cuzzi 
2006; Kennedy & Kenyon 2008). 

Another uncertain parameter is S, the disk surface den- 
sity power-law index. While we used S = 3/2, many 
models also consider S = 1. With S = 1, the main re- 
sults of Figure 2 are unchanged, with stronger migra- 
tion accounting for lower mass planets as the snow line 
distance decreases. It is unlikely observations of short- 
period super-Earths can constrain a or S based on Fig- 
ure 2, because they affect lower limits to planet masses, 
which will be hard to detect. 

The efficiency of type I migration is also unclear. Our 
choice of /„iig = 10 is based on numerical simulations, 
but may also be probed by future discoveries. The max- 
imum stellar mass M^^^max is not very sensitive to the 
snow line distance or disk profile, so for fixed snow line 
and metallicity enhancements (A), observations probe 
values for /mig. 

Our final major assumption is that planets form 
rapidly, and reach isolation before migrating. If plan- 
etesimals are small and growth is shear dominated, this 
assumption is generally true. With larger planetesimals 
however, growth is slower and planets may leave their 
formation regions while still growing. Planetary growth 
and migration are then coupled, and must be calculated 
simultaneously. Recently, Chambers (2006a, b) showed 
how a semi-analytic model of oligarchic growth can take 
different planetesimal sizes into account, and estimate 
their effect on growth rates (see also Thommes ct al. 
2003; Brunini & Benvenuto 2008). We now turn to a 
similar, yet simplified model to estimate the effects of 
planetesimal size on growth and migration. 

4.2. Semi-Analytic Model 



If planets grow fast enough, the isolation mass sets 
the range of disk masses that form migrating planets. If 
planetesimals are large enough, growth is not shear domi- 
nated and is slower. Migration then begins before planets 
reach isolation. To follow this evolution, a model treat- 
ing concurrent accretion and migration is necessary. Our 
model tracks damping of planetesimal random velocities 
by gas drag and stirring by a growing protoplanet. The 
random velocities set how growth proceeds relative to 
migration, allowing comparison with the analytic model. 

In the model, a single protoplanet of mass A/pi grows 
on a circular orbit from a planetesimal disk of small 
bodies of radius r. We adopt the accretion rate of 
Inaba et al. (2001) with the atmosphere enhanced ac- 
cretion radius of Inaba & Ikoma (2003). To account 
for accretion of other nearby protoplanets, the growth 
rate is increased by 50% (Chambers 2006b). Plan- 
etesimal random velocities are stirred by the growing 
protoplanet (Ohtsuki et al. 2002) and damped by gas 
drag (Inaba et al. 2001). The protoplanet accretes and 
stirs material within an annulus of half- width 4 Rn , 
and undergoes type I migration at the rate derived by 
Tanaka et al. (2002), modified by the offset /mig- We 
use a ten times less efficient migration rate, motivated 
by numerical (D'Angelo et al. 2002, 2003; Masset et al. 
2006a), and Monte-Carlo simulations (Ida & Lin 2008). 
Objects have mass density p = 1.5 g cm~"^ outside the 
snow line. Simulations are started with planetesimals 
in an equilibrium between protoplanet stirring and gas 
drag. Because we consider growth only near the snow line 
(see below), planetesimals do not undergo radial motions 
due to gas drag. Planetesimals lost to gas drag can be 
returned to the growth region by the cold finger mecha- 
nism (Cuzzi & Zahnle 2004). The system is evolved us- 
ing 4th order Runge-Kutta integration with an adaptive 
step-size (Press ct al. 1992). 

As before, we model protoplanets that form just out- 
side the snow line. These are the largest objects that 
migrate to the central star in a cold finger disk and are 
largely unaffected by interior objects. However, a mi- 
grating protoplanet shepherds material inward as it mi- 
grates, and will accrete some terrestrial material. This 
accretion cannot be treated by the semi-analytic model, 
so protoplanets cease accretion once they pass inside the 
snow line in the semi-analytic model. We model shep- 
herding with n-body simulations in §4.3. We vary j] to 
form 1-10 Af0 planets and use A = 10. 

Protoplanets begin with masses 1 x 10^^ M®, at 4i?H 
outside the snow line. This starting condition allows 
them to reach isolation if growth is faster than migra- 
tion. The disk is split into 1000 equally spaced radial 
bins. However, because accretion inside the snow line 
is turned off, objects grow from material in ~100 bins 
outside the snow line. 

The snow line distance and gas disk are as in §4.1 (eq. 
5 and following text), but the surface density of the gas 
disk decays exponentially with an e-folding time of 1 Myr. 
We place the inner edge of our disk at 0.2 M^, AU, though 
planets that reach a few tenths of an AU are migrating 
so rapidly that the exact value matters little. 

To test our code, we compare growth at 5 AU with 
Figure 1 from Chambers (2006b). His figure compares 
isolation times for different r with the type I migra- 
tion timescale. The smallest size planetesimals allow 
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Fig. 4. — Growth to isolation of a protoplanct at 5AU around 
a Solar-mass star. Isolation times arc in good agreement with 
Figure 1 of Chambers (2006b). Both models have cro = lOg cm~^ 
and tTg = 900 g cm ~^ . Our model includes explicit calculation 
of planetesimal eccentricities and inclinations, which accounts for 
differences. Lines are for r = 100 m, 1km, 10 km, and 100 km from 
left to right. The dashed (dotted) lines show the type I migration 
timescalc for /mig = 1 (10)- 



protoplanets to reach isolation before migration starts. 
Growth was simulated at 5 AU around a Solar-mass star, 
with a solid surface density of 10 g cm^ (so M\so ~ 
10 M0), and a gas/solids ratio of 90. Migration was not 
included, and the isolation time was simply compared to 
the analytic estimate of Equation (3). 

Figure 4 shows growth in the absence of migration at 
5 AU around a Solar-mass star for a range of r with 
similar initial conditions. The time to reach isolation is 
fastest for the smallest r (100 m), because growth is al- 
ways shear dominated. For r = 1 km, the growing proto- 
planet excites the small body random velocities. Growth 
ceases to be shear dominated at several 10^ yr. For higher 
r, isolation takes even longer, due to the decreasing effec- 
tiveness of gas drag on larger planetesimals. Compared 
with Figure 1 from Chambers (2006b) , the time to reach 
10 M^ is in good agreement. Our explicit calculation of 
eccentricities and inclinations accounts for differences in 
how growth proceeds (c.f. Fig. 3 of Chambers 2006a). 

Models with sufficiently small planetesimals reach iso- 
lation before migration. With /,„ig = 1, J' ^ 100 m, 
and for /,„ig = 10, r < 1 km. Thus, even with a reduced 
migration rate, protoplanets may still migrate before iso- 
lation. 

4.2.1. Semi- Analytic Model Results 

For the range of disk masses (ry) that forms 1-10 M^ 
planets around stars with masses 0.25-2 M0, Figure 5 
shows semi- major axis and mass evolution for r = 100 m. 
The choice of 1 M^ is somewhat arbitrary, but represents 
a rough lower limit for detection. We first describe the 
Solar case, and then look at differences as the stellar 
mass, and r change. 

For a Solar mass star, growth is not always fast enough 
for migration to occur before the gas disk is dispersed. 
For 77 = 2, the objects Hill radii increase faster than small 
bodies are stirred; thus growth remains shear dominated 
(^i* ^ !)■ The protoplanct successfully migrates to the 
inner edge of the disk. For lower 77, stirring overcomes 
damping at several xlO'^yr and growth slows. Higher 77 



results in faster growth of larger objects, which migrate 
early enough to avoid stalling at intermediate radii. With 
77 = 1, migration is somewhat significant, and the '^3 
Earth mass planet stalls at ~1 AU due to dissipation of 
the gaseous disk. For 77 = 0.5, the Earth-mass planet 
migrates little, and remains beyond the snow line. Final 
planet masses, and the degree of migration, are set by 
the initial surface density beyond the snow line. 

We turn now to trends across a range of stellar masses. 
Because of smaller snow line distances, migration is eas- 
iest for planets in the 1-10 Mq range around lower mass 
stars. Low-mass stars are the most likely to form these 
planets, because the range of disk masses that form them 
is much larger. For higher mass stars the more distant 
snow line makes migration unlikely for all but the most 
massive planets. Growth is driven out of the shear dom- 
inated regime more easily due to lower gas density at 
greater distances. This result confirms the maximum 
stellar mass M^^^max described above. As in the analytic 
model, M^.inax lies between 1-2 Mq with A = 10, be- 
cause no planet with a mass ^10 Mq migrates signifi- 
cantly for 2Mq. As stellar mass increases, the relative 
disk mass required to form 1-10 Mq planets decreases 
(eq.^ 7). 

Figure 6 shows how growth changes if planetesimals 
are larger. Models again have A = 10, but now the 
planetesimal radius r = 10 km. For larger planetesi- 
mals growth is easily stirred out of the shear dominated 
regime by the large objects for all stellar masses. The 
disk masses needed to reach the same range of planet 
masses are higher, because planets migrate out of the ac- 
cretion region before they reach isolation. For 0.25 Mq 
stars, the maximum 77 = 30 only just forms 10 Mq plan- 
ets. Again, Af^^max lies between 1-2 Mq, indicating that 
it is largely independent of planetesimal size. Though 
growth is slower, the results for r = 10 km are largely 
the same as 100 m, because the surface density can be 
increased to account for the slower growth. 

Figure 2 also includes results from the semi-analytic 
model, showing the range of planet masses that reach 
short-period orbits for a range of stellar masses. Models 
were run for A = 10, with r = lOm-lOkm and M^, = 
O.1-2M0. The upper limit decreases at 0.1 M© due to 
an upper limit on disk masses. Results from the analytic 
model are in good agreement. The difference in the lower 
limit arises because migration is faster at smaller radii, 
allowing smaller planets to reach the inner disk edge in 
the semi-analytic model. 

In summary, using a more detailed migration model 
yields results similar to the simple analytic treatment 
in §4.1. The inclusion of growth rates due to different 
planetesimal sizes adds another dimension due to differ- 
ent relative timescales for migration and accretion. The 
model offers more insight into how growth proceeds, and 
how the growth rate sets the required disk mass for form- 
ing short-period planets. 

4.3. Shepherding 

As a large body migrates inward, it captures interior 
objects onto mean motion resonances, and shepherds 
them inward. In the original scenario, a gas giant forms 
near the snow line, and subsequently migrates inward. 
As the giant migrates it shepherds interior protoplanets 
inward, which collide and merge to form super Earth- 
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Fig. 5. — Results from the semi-analytic cold finger disk model with A = 10 and r = 100m. Each vertical pair of panels shows mass 
and semi-major axis evolution for several relative disk masses. Higher disk masses (rj in legends) form more massive planets that migrate 
earlier. All 1— lOAfg planets migrate to the inner disk edge for 0.25 A^q, while none do for 2Mq. 



mass planets (Zhou ct al. 2005). Here we use n-body 
simulations to study a similar scenario, but with a low- 
mass planet migrating inwards from the snow line due to 
type I migration. 

To investigate shepherding effects, we used the MER- 
CURY integrator (Chambers 1999), including type I mi- 
gration and damping forces (Cresswell & Nelson 2006; 
Fogg & Nelson 2007). The migration rate, and eccen- 
tricity and inclination damping are reduced by a factor 
of /mig = 10. The inner disk edge is placed at 0.05 AU, 
and inside this point planets cease to interact with the 
disk (Terquem & Papaloizou 2007). Simulations are ini- 
tialised with a number of isolated protoplanets in a disk 
between 0.1 AU and the snow line at 2.7M^AU. 

Isolation masses are calculated from Equation (7) with 
the half-spacing B randomly varied between 3.75 and 
4.25. We assume Solar metallicity and /ice = 10. One 
protoplanet begins beyond the snow line. This outermost 
protoplanet is «30 times more massive than the one im- 
mediately interior to it (eq. 2). Initial eccentricities (in- 
clinations) are randomly distributed between and 0.02 
(0.5°), and the remaining orbital elements are randomly 
distributed. We set the mass of the outermost planet at 
the middle of the range shown for M^, = 0.25, 0.5, and 
1 Mq in Figure 2; 2, 3.2, and 6M0 respectively. Simula- 
tions are run for 10* yr with ^0.3 day timesteps. Objects 
are allowed to collide, and are assumed to merge into a 
single body with no fragmentation. These simulations do 
not include relativistic effects, or tidal interaction with 
the star. See Terquem & Papaloizou (2007) for a more 



detailed study of migration to small radii, and how these 
effects affect final system dynamics. 

4.3.1. Shepherding results 

Figure 7 shows the semi-major axis evolution result- 
ing from these simulations. All show similar charac- 
teristics. Starting from the inner disk edge, a wave 
of chaotic growth moves outward (e.g. Chambers 2001; 
Kenyon & Bromley 2006), until the number of proto- 
planets is reduced such that their spacing is stable. This 
stability is set by a balance between mutual perturba- 
tions between protoplanets, and damping by interaction 
with the gas disk. 

When the outermost large protoplanet begins to mi- 
grate, it scatters the first objects it encounters into ex- 
terior orbits. When the interaction occurs, these objects 
are still undergoing eccentric chaotic growth, and are less 
likely to be captured onto resonances and shepherded in- 
ward. Once scattered, the outer objects slowly migrate 
inward. For 0.25 and 0.5 M© the scattered planets are 
still relatively close to the star, and have time to set up 
chains of (mostly first order) resonant orbits. A few col- 
lisions occur. For IMq more objects are scattered out- 
ward, which continue chaotic growth. Despite the initial 
disruption by the migrating object, '^Earth-mass planets 
stiU form at ~1 AU. 

When the migrating protoplanet encounters objects 
that have reached stable orbits, it shepherds them in- 
ward. These smaller objects accrete others as their or- 
bits are pushed together, and several ~1 Mq rocky ob- 
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jects form. Shepherded objects may be accreted by the 
large migrator, or remain in interior resonant orbits (see 
Tcrqucm & Papaloizou 2007). 

These simulations show that as with the gas giant case 
(Zhou et al. 2005), the effect of super-Earth migration on 
interior objects has observational consequences. Planets 
near the outer edge of the terrestrial region are scattered 
outward, while those in the inner region are shepherded 
to smaller radii. Shepherding results in multiple short- 
period planets with different compositions. 

5. DISCUSSION AND SUMMARY 

We have considered two scenarios for forming short- 
period <1OM0 planets over a range of stellar masses: 
planet-planet scattering and type I migration. 

Our models form planets in cold finger disks. 
These disks have large snow line enhancements com- 
pared to the MMSN model (Stevenson & Luninc 1988; 
Cuzzi & Zahnle 2004). Water vapour from the terres- 
trial region condenses into ices outside the snow line as 
the gas disk diffuses and advects. The enhancement is 
increased by new water vapour delivered inside the snow 
line by drifting icy planetesimals (Cuzzi & Zahnle 2004). 
Protoplanets forming in the cold finger regions near the 
snow line are much larger than others elsewhere in the 
disk. 

We test the effectiveness of planet-planet scattering 
with n-body simulations. We consider stars with masses 
M^ = 0.25-2 Mq and 10 M^ planets. Planets with orbits 
near the limits of stability are evolved until a collision or 
ejection occurs, or 1 Gyr. Although equal mass planet- 



planet scattering can produce planets with small perias- 
tra for the lowest mass stars (Fig. 1), long circularisa- 
tion times prevent them from achieving circular orbits on 
reasonable timescales. Thus, scattering is probably not 
a viable scenario for placing low-mass planets on short- 
period orbits for any stellar mass. For 0.25 Mq, planets 
have periastra ~0.05AU and semi-major axes ~0.5AU. 
Though transit durations are still several hours, orbital 
periods of several hundred days and maximum radial ve- 
locities of a few m/s make these planets hard to detect. 

Migration of icy protoplanets from the snow line is a 
viable mechanism for forming short-period super-Earths. 
Planet masses set whether they migrate to the inner disk 
edge before the gas disk disperses. Some planets with in- 
sufficient masses arc stranded at intermediate radii as 
the gas disk disperses; a way to form "ocean planets" 
(Kuchner 2003; Leger et al. 2004). The minimum pro- 
toplanet mass for migration to a close-in orbit increases 
as the snow line moves out with increasing stellar mass 
(Fig. 2). The maximum planet mass is ^lOM^, because 
above this mass they instead accrete large atmospheres 
and form gas giants. 

Above ~ 1 Mq, the only protoplanets massive enough 
to migrate to close-in orbits are >10 Mg, and no hot 
super-Earths form. This maximum stellar mass is inde- 
pendent of the disk mass distribution, and probes type 
I migration efficiency. Other uncertain parameters, such 
as snow line distance and disk profile do not have major 
observable consequences, but are not easily constrained 
by observations either. 

For disks with similar snow line enhancements, the the- 
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ory yields trends with metallicity (Fig. 3). For disk 
masses distributed as a power law, the frequency of short- 
period planets increases with metallicity, because most 
disks have low masses. However, if disk masses arc dis- 
tributed around a relatively high mass, planet frequency 
decreases with increasing metallicity, because planets 
forming in the most common disks are pushed above the 
gas accretion mass at high metallicitics. As planetesimal 
size increases, growth slows, and becomes longer than the 
migration timescale. Simulations of concurrent accretion 
and migration with increased planetesimal sizes require 
much higher disk masses to yield similar results. 

As icy planets migrate from the snow line, they inter- 
act dynamically with interior rocky protoplancts (Fig, 
7). Protoplancts undergoing chaotic growth are scat- 
tered onto exterior orbits. Closer protoplancts on sta- 
ble orbits damped by disk interaction are shepherded in- 



ward, and coalesce into a few rocky objects with masses 
~1 M®. These objects may be accreted by the large mi- 
grating planet, or remain as separate planets on interior 
orbits. These orbits are likely ncar-commcnsurate with 
the icy migrators orbit (Tcrqucm & Papaloizou 2007). 
If planetary systems in such configurations are found in 
transit surveys, compositional models may discern dif- 
ferences, thus confirming their origins in rocky or icy 
regions (Valencia et al. 2007). However, different struc- 
tural models may be degenerate if the planets have at- 
mospheres (Adams et al. 2007). 

Some planets may accrete hydrogen atmospheres due 
to a decreased planetesimal accretion rate following iso- 
lation (Ikoma et al. 2000; Rafikov 2006). To be observed 
as hot super-Earths requires subsequent photoevapora- 
tion (e.g. Baraffe et al. 2005). Significant photoevapo- 
ration of planets with massive atmospheres is unlikely 
unless the planet mass is in the <7OM0 type I migra- 
tion regime (Raymond et al. 2007). Thus, planets with 
remnant hydrogen atmospheres may form by the same 
migration mechanism we present here. Scattering is also 
a possibility for these planets to reach short-period or- 
bits, because they have higher initial masses. 

For planets originating in icy regions, their largely 
volatile composition has important implications for their 
evolution during and after formation. Icy grains may 
enhance growth if they stick together more easily, but 
also allows the possibility of large evaporation events in 
high energy collisions of larger objects. During the vi- 
olent accretion process, and with the possible outcome 
of short-period orbits, melting and evaporation of ices 
will affect these planets (e.g. Lunine & Stevenson 1982; 
Kuchncr 2003; Selsis et al. 2007). 

After migrating to close-in orbits, initially icy /watery 
planets may retain large super-critical steam atmo- 
spheres, or become rocky cores stripped of volatiles 
entirely. Kuchncr (2003) considered the existence of 
volatile-rich planets in the Solar habitable zone, and sug- 
gested that planets around Solar luminosity stars would 
be safe from evaporation at >1 AU but not at closer dis- 
tances. He also noted that lower EUV luminosities for 
M dwarfs makes these stars less likely to evaporate plan- 
etary atmospheres. More recently, Selsis et al. (2007) 
revisited the issue, and concluded that planets >6Mj9 
will retain most of their water content at >0.04 AU from 
a Solar-type star. The results of both studies suggest 
the evaporation timescale is strongly dependent on semi- 
major axis. Therefore, a trend may be noticeable within 
the small semi-major axis range of transiting planets. 

The picture that emerges is of systems with evapo- 
rated rocky planets inside ~0.04AU, and steam plan- 
ets somewhat outside this distance. A few stalled ocean 
(Kuchncr 2003; Lcger et al. 2004) and icy planets extend 
through and past the habitable zone. For these planets, 
microlensing provides sensitivity complementary to tran- 
sit and radial velocity methods at ~AU distances (e.g. 
Beaulieu et al. 2006), which will help yield trends with 
semi-major axis, particularly for low-mass stars. 

Surveys such as the MEarth Project 
(Nutzman & Charbonneau 2007), CoRoT (Baghn 
2003), and Kepler (Borucki et al. 2003) hope to discover 
super-Earths by the transit method. Like those discov- 
ered by radial velocity, most planets will orbit close to 
their parent stars. Because they are unlikely to form 
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in situ, these planets necessarily require some form of 
migration or scattering from their formation regions. 
Observed systems will thus test and inform mechanisms 
that form and bring planets to visible orbits. 
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